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Neutron scattering experiments have revealed a field-induced antiferroquadrupolar order param- 
eter in the Pr-based heavy-fermion superconductor PrOs4Sbi2. We observed the field-induced anti- 
ferromagnetic dipole moment with the propagation vector q — (10 0) for the applied field direction 
both H |[ [110] and [0 01]. For JT||[110] at 8T, it should be noted that the induced antiferromag- 
netic moment of 0.16(10) (j,b/P? orients parallel to the field within our experimental accuracy. This 
observation is strong evidence that the O xy electric quadrupole to be the primary order parameter 
for JT||[110]. A mean-field calculation, based on Pi singlet and F^ 2 ' excited triplet with the O xy 
quadrupolar interaction, reproduces the induced moment direction and its field response. These 
facts indicate the predominance of O xy -type antiferroquadrupolar interaction in PrOs4Sbi2. 

PACS numbers: 75.40.Cx, 75.25.+Z, 61.12.Ld, 74.70. Tx, 75.30.Kz 



I. INTRODUCTION 

PrOs4Sbi2 is the first Pr-based heavy-fermion super- 
conductor and attracts considerable interest because of 
unusual superconducting properties ,ii2iMiSiSiLSiS One of 
the most interesting topics is what is the role of the 
multipole moments of Pr / electrons for the heavy- 
fermion nature and unconventional superconductivity in 
PrOs4Sbi2 i lc>l11 Long range magnetic dipole order does 
not coexist with superconductivity at zero fieldi^ In- 
stead, the existence of a magnetic field-induced ordered 
phase has been reported above the upper critical field H C 2 
for the superconductivit y 12 i 13 i 14 ' 15 In our previous study 
antiferromagnetic reflections with q=(10 0) were ob- 
served in the field- induced ordered phase for i? || [0 l]. 16 
The antiferromagnetic moment of 0.02 ^e/Pr at 8T 
along the [0 1 0] direction perpendicular to the applied 
field was significantly smaller than the field-induced fer- 
romagnetic moment parallel to H . This small dipole mo- 
ment and the H-T phase diagram can be well explained 
by the mean-field calculation based on the crystal field 

(2) 

level scheme of the Pi singlet ground state with T\ 
triplet excited state and assuming an antiferroquadrupo- 
lar (AFQ) interactionj 16 : 17 : 18 The field-induced ordered 
phase appears as a result of the level crossing of the zero 
field ground and excited state which comes from the field 
split triplet state. The level scheme and its field response 
have been experimentally clarified by means of neutron 
inelastic scattering under magnetic fields* 1 ^ The overall 
aspect of the field-induced ordered phase indicates that 
the antiferromagnetic dipole is not a primary order pa- 
rameter. 

Generally speaking, there are many possibilities for the 
primary multipole order parameters. The ordered phase 
is stable only under applied magnetic field, where the 
time reversal symmetry has been broken, thus, the anti- 



ferromagnetic dipole order appears simultaneously with 
the coupled order parameters such as quadrupole, oc- 
tupole, and hexadecapole in the order parameter space 
with given Th symmetry^ Since neutrons exhibit no 
cross section for an electric quadrupole moment, the ob- 
servation of the field-induced antiferromagnetic peak by 
neutron is not a direct proof for the primary AFQ order 
parameter. In this sense, our previous study did not rule 
out possibilities of other high-rank multipoleSf 1 ^ Actu- 
ally, it is shown that the qualitative features for £T|| [00 1], 
the induced antiferromagnetic moment perpendicular to 
the applied field and the phase diagram can be repro- 
duced in terms of octupole as well as quadrupol o 18 ' 20 . 
Recent intensive studies on skutterudites report an im- 
portant role of higher rank multipole, octupole^ 1 - and 
hexadecapole^ 2 , in some compounds. Therefore, it is in- 
dispensable to clarify strictly the dominant interaction 
in PrOs 4 Sbi2. For this purpose, further information is 
required: field direction and pressure dependencies and 
inelastic neutron scattering. 

Careful magnetization measurements on PrOs 4 Sbi2 
clarified that the field-induced ordered phase also exists 
for [1 1 0] and [1 1 1] with remarkable anisotropy in the 
H-T phase diagram^ The existence of the field-induced 
ordered phase for these different field directions also sup- 
ports the singlet-triplet crystal field level scheme. It 
has been pointed out that the anisotropy in the phase 
diagram is intrinsic to Th symmetry and the predomi- 
nance of quadrupolar interactionj 17 ' 18 The field-induced 
ordered states for both H\\[\ 1 0] and [1 1 1] are also con- 
sidered to be the AFQ ordered phase and not to be the 
dipolar one. The order parameter for H\\[l 1 0] was the- 
oretically predicted to be O xy ^, but has not been con- 
firmed so far. If this is the case, one could discriminate 
the primary order parameter of non-magnetic quadrupole 
from coupled magnetic octupole since the induced anti- 
ferromagnetic moment direction is expected to be dif- 
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FIG. 1: A photograph of the single crystalline sample of 
PrOs4Sbi2 used for the present neutron scattering study. 



ferent between them. Therefore the study for 1 0] 

is worthwhile in order to clarify a dominant interac- 
tion in PrOs4Sbi2- Besides, it motivates us that the 
magnetization measurements revealed the existence of 
some additional anomalies within the hcld-induced or- 
dered phase for 2T||[110], indicating possible transitions 
between multipolar ordered states»« 

The purpose of the present neutron scattering study 
is to clarify the order parameter of the field-induced or- 
dered phase for H || [1 1 0] and to investigate the origin of 
anomalies in this ordered state. We found that the in- 
duced antiferromagnetic moment is parallel to the field 
direction, which is the strong evidence for the O^-type 
quadrupole as the primary order parameter. 



II. EXPERIMENT 

Neutron scattering experiments have been carried out 
on the cold neutron triple-axis spectrometer LTAS in- 
stalled in the guide hall of the research reactor JRR-3 in 
Japan Atomic Energy Agency, JAEA. The instrumental 
setup was just the same as our previous experiment^ 
The use of a cold neutron with triple-axis mode was cru- 
cial for taking high quality data to reduce the contamina- 
tion of elastic spectra by the low-energy excitations. We 
used a liquid He- free 10 T superconducting magnet and 
3 He- 4 He dilution refrigerator both developed by JAEA. 23 
The vertical field was applied along [110], perpendicular 
to the (h h I) scattering plane. We observed superlatticc 
peaks at h + k + I — odd, while ferromagnetic scattering 
is superposed on nuclear Bragg peak at h + k + l = even, 
for example (110). 

A large single crystalline sample with a mass of 6 g 
has been grown by the antimony-self-flux method. The 
details of the sample preparation have been published 
elsewhere Figure Q] shows the picture of the sam- 
ple, which is composed of small crystallites typically 1 
mm 3 distributed within 1 degree of the crystallographic 
axes. The four-circle x-ray crystal structure analysis on 
the small piece of this sample confirms the filled skut- 
terudite structure with the space group Imi (T^, #204). 
The sample was characterized via magnetic susceptibility 
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FIG. 2: (a) The field dependence of the 2 21 superlattice re- 
flection profile of PrOs 4 Sbi2 measured at 0.29 K for H\\[l 1 0]. 
The solid lines denote Gaussian fitting, (b) The field de- 
pendence of the 2 21 reflection intensity measured at various 
temperatures, (c) The 1 1 reflection intensity plotted as a 
function of applied magnetic field taken at 0.27 K. The inset 
gives its differential curve, AI /AH. 



measurements; we observed clear Curie- Weiss behavior 
at high temperature and the maximum at 3.6 K, which is 
identical to the high-quality sample used in the de Haas 
van Alphen study£ The resistivity was also very similar 
to the previous results and our sample showed supercon- 
ducting transition at the onset of 1.88 K. These results 
guarantee the high quality of the present sample. For- 
tunately there was no domain structure concerning the 
two-fold symmetry in the basal plane, thus for example, 
(hkl) reflection is generally different from (khl). We do 
not know the mechanism for the single domain growth 
process, nevertheless the single domain sample can be re- 
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FIG. 3: (a) The temperature dependence of the 2 21 superlat- 
tice reflection of PrOs4Sbi2 measured under the magnetic field 
of 8 T applied along the [1 1 0] direction. The solid lines de- 
note the result of fitting with Gaussian and linear background 
for each temperature, (b) The temperature dependence of the 
2 21 reflection intensity measured at various fields. The solid 
lines are guides for the eyes. 



producibly obtained with our sample growth condition, 
which is rather mysterious. 



III. RESULTS 

Figure [U^a) shows the representative data for the scat- 
tering profile of the 2 21 superlattice peak measured at 
0.29 K with magnetic fields along the [110] direction. 
An application of the magnetic field of 8T induces a 
clear resolution-limited superlattice reflection at (2 21) 
where no trace of peak was observed at zero field. The 
observed field-induced superlattice peak positions are the 
same as those for H || [0 1], namely the propagation vec- 
tor is </=(100). A marked decrease to ^1/10 of the 
integrated intensity occurs between 8 and 10 T whereas 
the magnetic Bragg peaks maintain their positions and 
widths. 

The field dependence of the 2 21 peak intensity was 
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FIG. 4: The H-T phase diagram of PrOs 4 Sbi 2 for iJ||[l 10]. 
The closed and open circles are fields corresponding to the 
onset and the maximum of the field-induced antiferromag- 
netic peak, respectively. The closed triangles are the phase 
boundary determined from the temperature dependence. 



measured at various temperatures as shown in Fig. [21(b). 
Applying the field at 0.27 K, the peak intensity appears 
at 5T (=Ha), reaches a maximum around 8T (=-ff ma x), 
and decreases steeply above 8 T. Ha increases as the tem- 
perature rises: 6T and 7T at 0.51 K and 0.71 K, respec- 
tively. In contrast, the maximum around 8T is almost 
temperature independent. No trace of superlattice peak 
was found above 1 K. 

Figuredjc) displays a field variation of the peak inten- 
sity at 0.27K taken at Q=(l 1 0). An increase of 1 1 re- 
flection intensity originating from the uniform magnetic 
moment was clearly observed. The inset gives the differ- 
ential curve, dl/dH, which corresponds to the differen- 
tial curve of the square of the ferromagnetic moment M 2 . 
With applied fields, the intensity exhibits a gradual in- 
crease and changes its slope around 5.5 T and 8T, which 
can be clearly seen in the inset. These inflection fields 
are consistent with the onset (-Ha) and maximum fields 
(_ff max ) in the superlattice reflection intensity as well as 
the result of magnetization.—. The small hump around 
1.5 T might arise from the breaking of the heavy- fermion 
superconductivity since it is very close to the upper crit- 
ical field H C 2- 

The temperature dependence of the 2 21 superlattice 
peak profile under the applied field of 8 T is shown in 
Fig. [31(a). The 2 21 superlattice reflection becomes weak 
without significant peak shift and/or broadening as tem- 
perature increases, and no trace of peak was observed 
at 1.4 K. Figure [3](b) shows peak intensity at (2 21) as 
a function of temperature taken under the field of 3, 8 
and 10 T. With increasing temperature from 0.27 K at 
8 T, the peak intensity showed the monotonous decrease 
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FIG. 5: The product of the magnitude of the magnetic mo- 
ment and the angle factor derived from the observed magnetic 
Bragg reflection intensity. The open and closed circles indi- 
cate the data for hqH =8T and 10 T, respectively. The dot- 
ted lines are the fitting results for the reflection in the (hhl) 
scattering plane by assuming the antiferromagnetic moment 
parallel to the applied field, namely, sin0 = 1. 



and disappeared around Ta=0.8K. The feature for 10 T 
is quite similar to that for 8 T except for the absolute 
intensity. No peak was observed in the result for 3T. 
No additional anomaly at 8T and 10 T is consistent with 
the fact that -ff max is almost temperature independent. 
The present results for fT||[l 10] of PrOs4Sbi 2 are sum- 
marized in the H-T magnetic phase diagram as shown 
in Fig. 21 The closed and open circles are the field cor- 
responding to the onset and the maximum of the field- 
induced antiferromagnetic peak, respectively. The closed 
triangles are the phase boundary determined from the 
temperature dependence of the antiferromagnetic inten- 
sity. The overall aspect of the phase boundary, deter- 
mined from the present diffraction experiments, is con- 
sistent with the result of magnetization measurements^ 
The field-induced superlattice reflection with g=(10 0) 
exists only in the region of the field-induced ordered 
phase, in other words, the ordered state with q=(l 0) is 
the basis of the field-induced ordered phase for H ||[1 1 0] 
which is the same as fT||[0 01]. 

In order to determine the magnetic structure in the 
field-induced ordered phase, the integrated intensity of 
the superlattice peaks on the (hhl) scattering plane were 
measured. The magnetic Bragg reflection intensity for 
unpolarizcd neutron diffraction can be written, 

W<9) = AX(0)[O.269xl(T 12 /(Q)(/isiii0 | F mag (Q) |] : 

(1) 

where K is the scale factor, L(9) is the Lorentz factor, fx 
is the antiferromagnetic moment, f(Q) is the magnetic 
form factor and F ma _ g (Q) is the magnetic structure fac- 
tor which is unity for the present case. <p is the angle 




FIG. 6: The obtained magnetic structure for the field-induced 
ordered phase of PrOs4Sbi2 with application of the magnetic 
field along ffplO]. 



between the ordered magnetic moment and the scatter- 
ing vector Q. The square of the product of angle factor 
and antiferromagnetic moment /i can be given as follows; 

1 AiSm0 l= \j KL(0)[O.269 x 10-^f(Q)F m ,^QW 

(2) 

Figure [5] shows the square root of integrated intensity 
divided by scale factor, Lorentz factor, structure factor 
and magnetic form factor of the Pr 3+ free ion. Therefore 
the vertical axis in Fig. [5] is the product of the antiferro- 
magnetic moment fj, and the angle factor sin cj). As clearly 
seen in Fig. O this quantity is almost isotropic within the 
present scattering plane. Although the data obtained at 
10 T is much weaker, it is also isotropic. This isotropic 
nature of the left hand side of equation Eq. ([2]) indicates 
that the antiferromagnetic moment is perpendicular to 
the (hhl) scattering plane at both 8T and 10 T within 
our experimental accuracy. Thus the induced dipole mo- 
ment observed by neutron is parallel to the direction of 
the applied field along [110]. Using this angle factor, 
sin</>=l, the induced antiferromagnetic moment for 8T 
and 10T is deduced to be /iaf = 0.16(10) /j,b and 0.07 /is, 
respectively. The antiferromagnetic moment exhibits a 
strong reduction in the magnitude from 8T to 10 T. 

There is a large induced ferromagnetic moment parallel 
to the field, which is roughly ten times larger than that of 
the field- induced antiferromagnetic moment?^ The sum 
of the ferromagnetic and antiferromagnetic moment gives 
the magnetic structure in the field-induced ordered phase 
as shown in Fig. [5] The moment size is different for 
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the each magnetic sublattice due to the antiferromagnetic 
component parallel to the field direction; the magnetic 
moment at the corner of the unit cell is larger than that 
at the center, or vice versa. 



IV. MEAN-FIELD ANALYSIS 

We present a mean-field calculation for the present re- 
sults with H\\ [1 1 0]. In PrOs 4 Sbi 2 , the Ti singlet ground 
state with the T\ first excited triplet state at 8 K is well 
separated from the other excited state . 16 ' 25 ! 26 This means 
that the low-lying singlet-triplet levels are responsible fol- 
low temperature physical properties. Hence, we move on 
our discussions to this low-lying singlet-triplet subspace. 
Active multipoles for the singlet-triplet system with J=4 
in the Oh and Th symmetries are classified in ref.[l3 and 
Il8l According to the symmetry lowering from Oh to Th, 
the triplet states in Th are represented by linear com- 
binations of those in Oh , whereas the Ti singlet is not 
influenced. The mixing is parameterized by y, the coeffi- 
cient in the crystal field Hamiltonian reflecting the effect 
of Th, where y = corresponds to Oh^ 

The possible multipole order parameters under mag- 
netic fields which cause admixture are classified by sym- 
metry analysis. There are two possible symmetries for 
-H" || [1 10], Ti and T 2 . For Ti, the magnetic dipole J x + J y 
is mixed with J x — J y depending on the magnitude of 



y. The induced antiferromagnetic moment in the field- 
induced ordered phase was revealed to be J x + J y in the 
present experiment. This result indicates the order pa- 
rameter for H || [1 1 0] to be the T\ symmetry which is con- 
sistent with the theoretical prediction^ Furthermore, 
the induced moment parallel to applied fields within the 
experimental accuracy suggests the small y parameter 
which is also consistent with the previous results on neu- 
tron diffraction for Zf||[001] and a small anisotropy in 
the H-T phase diagram. 

Hereafter, we study to what extent the simple 
quadrupolar interaction model can reproduce the present 
observed results. Our analysis is based on the Hamilto- 
nian consisting of crystal field potential, Zeeman energy 
and AFQ interaction, 



H = H t 



CF 



Hz + He 



(3) 



It is shown that the mechanism of field-induced AFQ 
order in this model can be clarified with a pseudo-spin 
representation analogous to the dimer-spin systems! 17 ' 18 
With the use of two pseudo-spins, s\ and s 2 with S — 
1/2, which are defined in the singlet-triplet subspace, the 
local part of the Hamiltonian is described as 

H C f + H z = I \ + s u ■ s 2l 

i ^ 

- h -^2(s u + s 2i + 5{s u - s 2i)), (4) 
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FIG. 7: The field dependence of the field-induced ferromag- 
netic and antiferromagnetic moment. The open circles are ob- 
tained from the neutron diffraction, whose value corresponds 
to the left axes. The results of mean-field calculation are rep- 
resented by lines corresponding to the right axes. 



where A is the singlet-triplet energy splitting and h is 
a scaled magnetic field. The appearance of an effective 
staggered field (5) to the two spins is characteristic of 
the Th system. It has been shown that the interaction 
of Ozy-type quadrupoles is mapped to the biquadratic 
pseudo-spin Hamiltonian, 

H Q = D Q 5Z 4 [( S u x 8 2i ) ■ ( s y x s 2 j) + £l/i i ' V i 

(ij) 

+ t2\ (s U X S 2i) ■ 11 



where 



s 1 s 2 , s 1 s 2 



Si s 



1.; 



1°2J °1°2 



2j. 



(5) 



(G) 



See ref. [IT] and [TH for details on the dependence of ei 
and e 2 on the Th parameter y. An advantage of this 
mapping is to clarify that the predominant and isotropic 
contribution to the field-induced order is given by the 
part of a vector product s i x s 2 , whereas the remaining 
part related with fi produces an anisotropic correction 
due to the Th symmetry. 

In the following, we discuss the mean-field solutions 
of this model in the field direction H\\[l 10], in which 
the stable phase is shown to be of O xy . The calculation 
has been carried out by using the same parameter set as 
in ref. 18. Namely, we assume Dqz — 0.3A, A = 8K, 
x = 0.45 and y = 0.12, where z is the number of the near- 
est neighbor bonds, and (a;, y) are the crystal field param- 
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FIG. 8: Calculated angle £ between the field-induced antifer- 
romagnetic moment and the field direction [1 1 0] as a function 
of relative magnetic interaction strength r\ for /i/A=1.2, cor- 
responding to /lo-ff— 7T. 

eters. The results of magnetic field dependence of the fer- 
romagnetic (Mp) and antifcrromagnetic moment (Maf) 
are shown in Fig. [7] They are compared with the present 
experimental results that are deduced by normalizing Mf 
and Maf at 8 T to be and 0.16/^b, respectively. 

One finds that the observed overall field response of both 
ferromagnetic and antiferromagnetic moments is well ex- 
plained by this model analysis. These agreement clearly 
indicates that the staggered moment should be induced 
by the primary AFQ order which remains inaccessible to 
the neutron probe. It is also shown that the calculated 
value of antifcrromagnetic moment is larger than that 
for if 1 1 [001], and this tendency is qualitatively consis- 
tent with the observed results. Note, however, that there 
exist quantitative differences between theory and experi- 
ment, probably due to fluctuation effects or contribution 
of other multipole interactions. 

In order to clarify the unique properties of the AFQ 
model more definitely, let us study the effects of pos- 
sible magnetic interaction. It is shown that the mag- 
netic interaction of dipoles and octupoles is expressed as 
an isotropic bilinear spin HamiltonianJ^ Although the 
generic form of the interaction is somewhat complex with 
a few parameters, we introduce here the simplest form to 
capture the essential physics, 

H M = D M ^2(s u ■ s ij + s 2i ■ s (7) 

This symmetric model in the two spins corresponds to a 
mixed interaction of dipoles and octupoles. Note, how- 
ever, that the relevant part is the octupolar interaction, 
because they have much larger matrix elements between 
singlet and triplet states. It should be stressed that the 
total Hamiltonian in the limit, Dq — and Dm ^ 0, re- 
duces to the conventional dimer-spin models, in which 
a similar field-induced order takes place by the magnetic 



interaction. Thus, it is quite instructive to see the differ- 
ence between quadrupolar and magnetic interactions by 
interpolating the two limits. 

In the framework of the mean-field theory, we have cal- 
culated the antiferromagnetic moment direction for the 
total Hamiltonian H+Hm in the field direction H\\[l 10]. 
Taking the ratio of the interactions as Dq — D(l — rf) 
and Dm — Drj, the angle from the applied field direction, 
C, is plotted in Fig. [5] as a function of relative magnetic 
strength r\ — Dm/(Dq + Dm) with a fixed Dz(— 0.3A). 
In the case that the magnetic interaction is small, the 
moment direction is almost parallel to the field. Whereas 
the magnetic interaction becomes dominant, the moment 
orients perpendicular to the field, as expected in the 
dimer-spin systems. Such a remarkable change of the 
antifcrromagnetic structure that depends on the inter- 
action is characteristic in the field direction if || [110]. 
In the case of if||[001], both interactions result in the 
antiferromagnetic moment perpendicular to the field. 30 
Thus, the present experiment for if || [110] is crucial in 
distinguishing the type of multipolar interactions. In this 
sense, the observed antiferromagnetic moment parallel to 
the field is regarded as strong evidence of the predomi- 
nant AFQ interaction. 



V. DISCUSSION 

The field-induced antiferromagnetic moment parallel 
to the field for if || [110] is different from the perpen- 
dicular antiferromagnetic moment for if || [0 1], whereas 
the propagation vector q=(10 0) is the same. In the lat- 
ter case the antiferromagnetic moment is parallel to the 
[0 1 0] direction, and not to the [1 0] direction; a clear 
peak was observed for q=(l 0) but no trace of the anti- 
ferromagntic peak was observed for q—(0 1 0). This char- 
acteristic moment direction reflecting the lack of four-fold 
symmetry gives the important information on the possi- 
ble order parameter. Variation in the orientation of the 
field-induced antiferromagnetic moment according to the 
applied field direction is characteristic to the quadrupolar 
ordered phase. Both results indicate the predominance 
of the O xy -type AFQ interaction in PrOs4Sbi2- 

A recent inelastic neutron scattering study clarified 
that the crystal field excitations from the Ti ground state 

(2) 

to the triplet first excited state have a dispersion; 
the observed excitation softens at q — (10 0) which is 
the same as the ordering vector of the field-induced AFQ 
ordered stated In addition, the intensity at q=(10 0) 
is weaker than that at zone center. Though the mag- 
netic and non-magnetic interaction gives the same en- 
ergy dispersion, the q dependences of the intensity are 
opposite^ namely, magnetic interaction should lead to 
the stronger intensity at the zone boundary. Thus, the 
observed excitation spectra indicate the dominance of an- 
tiferroquadrupolar interaction which is quite consistent 
with the present result for if ||[1 10]. Together with the 
anisotropy in the H-T phase diagram, the predominance 
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of AFQ interaction is clarified in PrOs4Sbi2- 

Furthermore, the excitation observed in the inelastic 
scattering study gives an interaction strength. Note that 
the quadrupolar coupling constant obtained from the in- 
elastic scattering is almost consistent with that used in 
the present study. There is a slight difference in the def- 
inition; Do in the present paper can be connected with 
c?q in refl25l with Dq = (3 2 <1q where (3 is the off-diagonal 
matrix elements of quadrupolar moments between sin- 
glet and triplet. Taking this difference into account, the 
coupling constant /3 2 dqz~-2.7K derived from inelastic 
spectra is consistent with Dqz~2.4K used in the present 
study. In addition, the quadrupolar coupling constant 
<?P 4 obtained from the ultrasound measurements of C44 
is also consistent with the present value.— In other words, 
the coupling constant used in the present study corre- 
sponds with those obtained in the other measurements. 

The maximum in the induced antiferromagnetic mo- 
ment, -ffmaxi corresponds well to the Hi anomaly in the 
magnetization^, as shown in Fig.0J However, no change 
in the antiferromagnetic structure was found on pass- 
ing through fl mal . As shown in Fig. the induced 
Jx + Jy appears to vanish slightly above 10 T. This crit- 
ical field is well below the reported critical field for the 
field-induced ordered phase of ~12T but rather close to 
the H2 anomaly^ In other words, there is another field- 
induced ordered state above Hi- A possible explanation 
would be the switching of the order parameter from T± to 
T2 representation at H2 . It is theoretically predicted that 
the energy difference of these two ordered states is very 
small. 18 The coupled magnetic dipole order parameter of 
T2 representation is J z for H\\[l 10]. Thus a change of 



moment direction from [1 10]||iJ to [0 01]_Li/ should be 
accompanied by the quadrupole-quadrupole transition. 
This reorientation of the antiferomagnetic moment has 
not been observed in our present neutron scattering ex- 
periment up to 10 T. In order to clarify the origin of the 
anomalies observed in magnetization, higher field neu- 
tron scattering experiments should be carried out. 

In conclusion, we observed the field-induced anti- 
ferromagnetic order of heavy-fermion superconductor 
PrOs4Sbi2 for H ||[1 1 0]. The antiferromagnetic moment 
was parallel to the field direction, which is direct evi- 
dence for the primary quadrupolar order parameter of the 
field-induced ordered phase. Together with our previous 
study for iT||[001], we can conclude that the O xy -type 
quadrupolar interaction is intrinsic to the many body 
interaction between / electron based quasiparticles in 
heavy-fermion superconductor PrOs4Sbi2- 
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In the case of the Oh symmetry (y — 0), the dipoles 
and the octupoles in the singlet-triplet system are sim- 
ply expressed with the pseudo spins as J as 1 + S2 and 
T |3 asi-S2.- When the T h term is introduced (y=£0), 
they are weakly mixed with each other. Concerning the 
mean-field solution of the interaction (7), it is known that 
the antiferromangetic order parameter is transverse to the 
field direction^. If the y component of the spins, or equiva- 
lently that of the dipoles and the octupoles, is chosen as the 
order parameter in H || [0 1], one can expect (s^) ~ — (sjj) 
due to the positive A term in eq. (4) . This leads to a tiny 
antiferrmagnetism because J y ~ s\ + s\ for a small Th 
term, which is compatible to the experimental results for 
H 1 1 [0 1]. Since the model is isotropic for the field direc- 
tions, this transverse antiferromagnetism does not account 
for the present experimental results in H ||[1 1 0]. 



